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' Consider the high-order heat-type equation du/dt = ±d N u/dx N for an integer 

N > 2 and introduce the related Markov pseudo-process (X(t))t^o- In this paper, we 
study several functionals related to (X(t)) t ^o: the maximum M(t) and minimum m(t) 
up to time t; the hitting times t+ and r~ of the half lines (a, +00) and (—00, a) respec- 
tively. We provide explicit expressions for the distributions of the vectors {X(t) 1 M(t)) 
and (X(t),m(t)), as well as those of the vectors (t+, AT(t+)) and (t~ ,X(t~)). 
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1 Introduction 

Let N be an integer greater than 2 and consider the high-order heat-type equation 

du d N u 

dt =KN d^ (L1) 

where kn = (— 1) 1+N / 2 if N is even and kn = ±1 if iV is odd. Let p(t;z) be the 
fundamental solution of Eq. (1.1) and put 

P(t;x,y) =p{t;x-y). 

The function p is characterized by its Fourier transform 

i 

e iu ^p(t;0dC = e KNt{ - iu)N . (1.2) 



/+oo 
-oo 



With Eq. (1.1) one associates a Markov pseudo-process (X(t))t^o defined on the real line 
and governed by a signed measure P, which is not a probability measure, according to the 
usual rules of ordinary stochastic processes: 

¥ x {X(t) G dy}=p(t;x,y)dy 
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and for = to < t\ < ■ ■ ■ < t n , xq = x, 

n 

F x {X(ti) G dxi, . . . ,X(t n ) G dx n } = - - Xi) dxi. 

Relation (1.2) reads, by means of the expectation associated with P, 




: iux+KNt(iu) 



Such pseudo-processes have been considered by several authors, especially in the partic- 
ulary cases N = 3 and N = 4. The case N = 4 is related to the biharmonic operator 
d A /dx A . Few results are known in the case N > 4. Let us mention that for N = 2, the 
pseudo-process considered here is a genuine stochastic process (i.e., driven by a genuine 
probability measure), this is the most well-known Brownian motion. 

The following problems have been tackled: 



• Analytical study of the sample paths of that pseudo-process: Hochberg [8] defined 
a stochastic integral (see also Motoo [14] in higher dimension) and proposed an 
Ito formula based on the correspondence dx A = dt, he obtained a formula for the 
distribution of the maximum over [0, t] in the case N = 4 with an extension to the 
even-order case. Noteworthy, the sample paths do not seem to be continuous in the 
case iV = 4; 

• Study of the sojourn time spent on the positive half-line up to time t, T(t) = 
meas{s G [0,i] : X(s) > 0} = J * l{x(s)>o} <^ s: Krylov [11], Orsingher [20], Hochberg 
and Orsingher [9], Nikitin and Orsingher [16], Lachal [12] explicitly obtained the dis- 
tribution of T(t) (with possible conditioning on the events {X(t) > (or =, or <)0}). 
Sojourn time is useful for defining local times related to the pseudo-process X, see 
Beghin and Orsingher [1]; 

• Study of the maximum and the minimum functionals 

M(t) = max X(s) and m(t) = min X(s) : 

Hochberg [8], Beghin et al. [2, 3], Lachal [12] explicitly derived the distribution of 
M(t) and that of m(t) (with possible conditioning on some values of X(t)); 

• Study of the couple (X(t), M(t)): Beghin et al. [20] wrote out several formulas for 
the joint distribution of X(t) and M{t) in the cases N = 2> and iV = 4; 

• Study of the first time the pseudo-process (X(t))t^o overshoots the level a > 0, 

= inf{t ^ : X(t) > a}: Nishioka [17, 18], Nakajima and Sato [15] adopt a dis- 
tributional approach (in the sense of Schwartz distributions) and explicitly obtained 
the joint distribution of r+ and X(t£) (with possible drift) in the case TV = 4. The 
quantity X(t^) is the first hitting place of the half-line [a, +oo). Nishioka [19] then 
studied killing, reflecting and absorbing pseudo-processes; 

• Study of the last time before becoming definitively negative up to time t, 0{t) = 
sup{s G [0, t] : X(s) > 0}: Lachal [12] derived the distribution of 0(t); 

• Study of Equation (1.1) in the case N = 4 under other points of view: Funaki [6], and 
next Hochberg and Orsingher [10] exhibited relationships with compound processes, 
namely iterated Brownian motion, Benachour et al. [4] provided other probabilistic 
interpretations. See also the references therein. 
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This aim of this paper is to study the problem of the first times straddling a fixed level 
a (or the first hitting times of the half-lines (a, +00) and (—00, a)): 

T + = inf{i ^ : X(t) > a}, r~ = inf{t ^ : X(t) < a] 

with the convention inf(0) = +00. In the spirit of the method developed by Nishioka 
in the case N = 4, we explicitly compute the joint "signed-distributions" (we simply 
shall call "distributions" throughout the paper for short) of the vectors (X(t), M(t)) and 
(X(t),m(t)) from which we deduce those of the vectors (r+, X(r+)) and (r~, X(t~)). The 
method consists of several steps: 

• Defining a step-process by sampling the pseudo-process (X(t))t^o on dyadic times 

t n>k = k/2 n , G N; 

• Observing that the classical Spitzer identity holds for any signed measure, provided 
the total mass equals one, and then using this identity for deriving the distribution 
of (X(t ni k), maxo^-^A; X(t n j)) through its Laplace-Fourier transform by means of 
that of X(t n! k) + where x + = max(x, 0); 

• Expressing time r+ (for instance) related to the sampled process (X(t n! k))k&N by 
means of (X(t nj k), max ^^ X(t n j)); 

• Passing to the limit when n — ► +00. 

Meaningfully, we have obtained that the distributions of the hitting places X(t£) and 
X{t~) are linear combinations of the successive derivatives of the Dirac distribution S a . 
In the case N = 4, Nishioka [17] already found a linear combination of 5 a and 5' a and 
called each corresponding part "monopole" and "dipole" respectively, considering that an 
electric dipole having two opposite charges S a + £ and 5 a - £ with a distance e tending to 
may be viewed as one monopole with charge S' a . In the general case, we shall speak of 
"multipoles" . 

Nishioka [18] used precise estimates for carrying out the rigorous analysis of the pseudo- 
process corresponding to the case N = 4. The most important fact for providing such 
estimates is that the integral of the density p is absolutely convergent. Actually, this 
fact holds for any even integer N. When N is an odd integer, the integral of p is not 
absolutely convergent and then similar estimates may not be obtained; this makes the 
study of X very much harder in this case. Nevertheless, we have found, formally at least, 
remarkable formulas which agree with those of Beghin et al. [2, 3] in the case N = 3. 
They obtained them by using a Feynman-Kac approach and solving differential equations. 
We also mention some similar differential equations for any N. So, we guess our formulas 
should hold for any odd integer N ^ 3. Perhaps a distributional definition (in the sense 
of Schwartz distributions since the heat-kernel is locally integrable) of the pseudo-process 
X might provide a properly justification to comfirm our results. We shall not tackle this 
question here. 

The paper is organized as follows: in Section 2, we write down general notations 
and recall some known results. In Section 3, we construct the step-process deduced from 
(X(t))t^o by sampling this latter on dyadic times. Section 4 is devoted to the distributions 
of the vectors (X(t),M(t)) and (X(t),m(t)) with the aid of Spitzer identity. Section 5 
deals with the distributions of the vectors (t+,X(t+)) and (t~ ,X(t~)) which can be 
expressed by means of those of (X(t),M(t)) and (X(t),m(t)). Each section is completed 
by an illustration of the displayed results therein to the particular cases N G {2,3,4}. 

We finally mention that the most important results have been announced, without 
details, in a short Note [13]. 
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2 Settings 



The relation f_°?p(t;£)d£ = 1 holds for all integers N. Moreover, if N is even, the 
integral is absolutely convergent (see [12]) and we put 

/+oo 
\p(t;0\dti>l. 
-oo 

Notice that p does not depend on t since p(t;£) = r 1 / iV p(l ; ^/t 1 /^). For odd integer N, 
the integral of p is not absolutely convergent; in this case p = +oo. 



2.1 N th roots of k n 

We shall have to consider the N th roots of kn (Oi for ^ / ^ N — 1 say) and distinguish 
the indices I such that 3?6>; < and > (one never has $l6i = 0). So, let us introduce 
the following set of indices 

J = {/ G {0,... ,JV- 1} : Wi > 0}, 
K = {/ G {0,... ,N- 1} : ^ < 0}. 

We clearly have J U K = {0, . . . , TV - 1}, J n if = and 

#J + #K = N. (2.1) 

If AT = 2p, then k w = (-l) p+ \ &i = e ^ 2l+p+l ^l N \ 

J = {p, ...,2p- 1} and K = {0, . . . ,p - 1}. 
The numbers of elements of the sets J and if are 

#J = #K = p. 

If iV = 2p + 1, two cases must be considered: 

• For k n = +1: 0/ = e^ 2 ^! and 

J={0,...,|}u{^ + l,...,2p} and K={| + 1,...,^} ifpiseven, 

J={0,..,^}u{^,..,2p}andK={^,...,^±l} ifpisodd. 

The numbers of elements of the sets J and K are 

# J = p + 1 and #if = p if p is even, 

# J = p and #if = p + 1 if p is odd; 

• For kn = -1: 0j = e i [( 2 ' +1 ) 7r / JV ] and 

J = {0, . . . , | - l} U + 1, . . . , 2p} and K = {|, . . . , ^} if p is even, 

J={0,...,^}u{^±l,...,2p}andK={^±i,...,^} ifpisodd. 

The numbers of elements of the sets J and if are 

# J = p and #if = p + 1 if p is even, 

# J = p + 1 and j^K = p if p is odd. 



Figure 1 illustrates the different cases. 
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J : Wj > K : W k < 

#2p-l 




J : »0,- >0 K : K0 fe < 



W2 




J : Mj > 



*(p-l)/2 



^3p/^ 3p/2+1 t7 (3 P +l)/2 ^(3p+3)/2 

Case iV = 2p + 1, kjv = +1: even p (left), odd p (right) 



K : < 

Op/2 



J : W-j > X : m < 




J : sftflj > 



(P"l)/2 




?3 P /2 *3 P /2+l (3p - 1)/2 l^(3p+l)/2 

Case N = 2p + 1, kn = — 1: even p (left), odd p (right) 
Figure 1: The N th roots of re^r 

2.2 Recalling some known results 

We recall from [12] the expressions of the kernel p(t; £) 



-i(u+K N t(-iu) N ^ u 



(2.2) 



together with its Laplace transform (the so-called A-potential of the pseudo-process (AT(t)) t ^o), 
for A > 0, 



r+oo 

$(A;0= / e- M p(t;0dt=< 
Jo 



1 

"iV 



A 1 /"" 1 ^^* for^O, 



i. A i/iV-i^^. e ^ for^O. 
je J 



keK 



(2.3) 



Notice that 



f+OO 

<&(A; = / e~ A * dt F{X(t) € 
JO 

We also recall (see the proof of Proposition 4 of [12]): 
*(A; = / e" A * P{X(t) < -f } rft = < 

JO 



k&K 



(2.4) 



for £ ^ 0. 



We recall the expressions of the distributions of M(t) and m(i) below. 
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Concerning the densities: 

c+oo 



with 



and 



/ e' u dtF x {M(t) G dz}/dz = - ip x {x - z) for x < z, 
Jo * 

/ e~ xt dt¥ x {m(t) G dz}/dz = -tp x (x-z) forx^z, 
Jo * 

va(o = E m,- V>a(£) = - E e ^ 

jGJ keK 



(2.5) 



(2.6) 



^= n 



9i 



forjGJ, 5 fe = -A- forfeGK 



6,1 " 



leK\{k} 



01 - 0k 



• Concerning the distribution functions: 



f 

Jo 



e~ xt ¥ x (M(t) < z)cft 



+oo 



-Ai 



; (m(i) ^ z) eft 



l-£B fce «*^(*-*) 



for x ^ z, 
for x ^ z. 



(2.7) 



We explicitly write out the settings in the particular cases N G {2,3,4} (see Fig. 2). 



Example 2.1 Case N = 2: we have k 2 = +1, 0o = -1,6*1 = 1, J = {l},if = {0}, 
Ai = l,B = l- ■ 



Example 2.2 Case iV = 3: we split this (odd) case into two subcases: 



• for n 3 = +1, we have 6 = l,0i = e i2w /_ 3 ,6 2 = e"* 2 ^ 3 , J 
Ao = 1,B 1 = TZ ^ ri =^e-^,B 2 = B 1 = ^ e^; 



{0},K = {1,2}, 

for k 3 = -1, we have 6 = e^l\O x = -1,6 2 = e~ i7T / 3 , J = {0,2}, K = {1}, 
Ao = tz^73 = 73 A 2 = A = -L e~^l\ B 1 = 1. 



Example 2.3 Case iV = 4: we have k a = -1, 6» = e i37r / 4 ,0i = e^ 3 ^ 4 , <9 2 = e~ in / 4 ,e 3 = 
e^/ 4 , J = {2, 3}, K = {0, 1}, A 2 = B = = 4- e~^'\ A 3 = B l = A 2 = j. e^/ 4 . 



2.3 Some elementary properties 

Let us mention some elementary properties: the relation n^i^l — e l ^ 2ln ^ N ^) = N entails 



N-l 

n 

=0,l^m 



0i _J_ 

0m - 01 N 



for < m € N - 1. 



(2.8) 
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N = 2 N = 3, k 3 = +1 JV = 3, k 3 = -1 N = 4 



Figure 2: The N th roots of k n in the cases TV G {1, 2, 3} 



The following result will be used further: expanding into partial fractions yields, for any 
polynomial P of degree degP ^ #J, 



P(z) 



^p(^) 



if degP<#J-l, 



A j p ( 9 j) + TT 5/ if deg P = # J and the highest 

gree coefficient of P is 1. 



(2.9) 



• Applying (2.9) to x = and P = 1 gives A? = YlktK^k = 1- Actually, the 

Aj's and P^'s are solutions of a Vandermonde system (see [12]). 



Applying (2.9) to x = 9k, k G if, and P = 1 gives 



E ft; A? 



9k/ 9 j 



X\(^-0k/9 j ) 



n -l 



N-l 

n 

n Si— fc 

;=0,;^fc 
iSK\{fc} 



which simplifies, by (2.8), into (and also for the Pfc's) 



0, - ft 



9j — 9k NP>k 



for k G K and — ^~ 

fcG_ftT 



(9 fc - 6>j JVA, 



for j G J. (2.10) 



Applying (2.9) to P = x p , p ^ #J, gives, by observing that l/0j = 9j, 



x>' 



\[{l-9 jX ) 



ifp<#J-l, 



— ^^ + (-i)# j - i rr^ ifp = #j. 



(2.11) 



3 Step-process 



In this part, we proceed to sampling the pseudo-process X = (X(t))t^o on the dyadic times 
t n ,k = k/2 n , k, n G N and we introduce the corresponding step-process X n = (X n (t)) t ^o 
defined for any n G N by 

oo 
k=0 
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The quantity X n is a function of discrete observations of X at times t njk , k G N. 

For the convenience of the reader, we recall the definitions of tame functions, functions 
of discrete observations, and admissible functions introduced by Nishioka [18] in the case 
N = 4. 

Definition 3.1 Fix n G N. A tame function is a function of a finite number of ob- 
servations of the pseudo-process X at times t n j, 1 ^ j ^ k, that is a quantity of the 
form F nk = F(X(t nj i), . . . , X(t rhk )) for a certain k and a certain bounded Borel function 
F : R fc — > C. The "expectation" of F n ^ k is defined as 



E 



x(F n ,k) = /••• / F(x 1 ,...,x k )p(l/2 n ;x,x 1 )---p(l/2 n ;x k -i,x k )dx 1 ---dx k . 

J JR k 



We plainly have the inequality 

\®x(F n , k )\ ^/sup|F|. 

Definition 3.2 Fix n G N. A function of the discrete observations of X at times t U:k , 
k ^ 1, is a convergent series of tame functions: Fx„ = YlkLi F n ,k where F n<k is a tame 
function for all k ^ 1. Assuming the series YlkLi \^x(Fn,k)\ convergent, the "expectation" 
of Fx„ is defined as 

oo 

^ x (F x J = Y,^x(F n , k ). 
k=i 

The definition of the expectation is consistent in the sense that it does not depend 
on the representation of Fx„ as a series (see [18]): if Ylk^=i Fn,k = Yl'kLi G n ,k and if 
the series 52T=i \ E *( F n,k)\ and 'EkLi\ E x{G ntk )\ are convergent, then Efeli E ^( i? n,fe) = 

^2T=i K x(G ntk ). 

Definition 3.3 An admissible function is a functional Fx of the pseudo-process X which 
is the limit of a sequence (F Xn )neN of functions of discrete observations of X: 

F x = lim F Xn , 

n^oo 

such that the sequence (K x (Fx n ))n£N is convergent. The "expectation" of Fx is defined as 

E X (F X ) = lim E x (F Xn ). 

n— >oo 

This definition eludes the difficulty due to the lack of a-additivity of the signed measure P. 
On the other hand, any bounded Borel function of a finite number of observations of X 
at any times (not necessarily dyadic) ti < ■ ■ ■ < t k is admissible and it can be seen that, 
according to Definitions 3.1, 3.2 and 3.3, 

E x [F(X(t 1 ),...,X(t k ))] = /•••/ F(x 1 ,...,x k )p(t 1 ;x,x 1 )p(t 2 -t 1 ;xi,x 2 )--- 

J JR k 

x p(t k - t k -i; x k -!,x k ) dx 1 --- dx k 

as expected in the usual sense. 

For any pseudo-process Z = (Z(t)) t ^o, consider the functional defined for A G C such 
that sft(A) > 0, (i G R, v > by 

f+oo 

F z (\,fi,v)= e - xt+i " Hz ^- vKz ^I z (t)dt (3.1) 
Jo 



where Hz, Kz, Iz are functional of Z denned on [0, +00), Kz being non negative and Iz 
bounded; we suppose that, for all t ^ 0, Hz(t), Kz(t), Iz(t) are functions of the continuous 
observations Z(s), ^ s ^ t (that is the observations of Z up to time t). For Z = X 
have 



wc 



OO „; 

Fx n (A, v) = Y, 



e -\t+i^H Xn {t n , h )- V K Xn (£ n fe ) dt 

00 / ft n fc+1 \ 

= Y( e-*dt\e i ' lHx »( tn ^- vKx »( tn ^Ix n (t ntk ) 

k=0 K 

1 - e"^ 



jfc=0 

00 / »f 1 1 1 



A 



. e -^t n ,k+i^H Xn (t nik )-vK Xn (t n ,k) i x (t nk ). (3.2) 



fc=0 



Since H Xn (t n:k ), K Xn {t n>k ), Ix n (t n ,k) are functions of X n (t nJ ) = X(t nJ ), ^ j ^ fc, the 
quantity e l ^ Hxn ^ n ' k ^~ uKx ^ tn ' k ^ Ix n (t n ,k) 1S a tame function and the series in (3.2) is a 
function of discrete observations of X. If the series 

00 

fc=0 

converges, the expectation of Fx n (A, /x, v) is defined, according to Definition 3.2, as 

^ E x \eT Xtn • k +W Hx n (*n,fc ) ~ ^X n (*n, fc ) J-y 



E^FxJA,/^)] 



fc=0 



Finally, if lim n ^ +00 Fx n (A, /U, z^) = Fx (A, /x, 1/) and if the limit of E x [F Xn (A, fi, v)\ exists as 
n goes to 00, -Fx (A, /x, zv) is an admissible function and its expectation is defined, according 
to Definition 3.3, as 

E x [F x (X,^,u)] = lim E x [F Xn (X,fi,v)]. 

rj._i._L.fv-i 



4 Distributions of p_(*),M(*)) and (X(t),m(t)) 

We assume that A" is even. In this section, we derive the Laplace-Fourier transforms of 
the vectors (X(t),M(t)) and (X(t),m(t)) by using Spitzer identity (Subsection 4.1), from 
which we deduce the densities of these vectors by successively inverting — three times — the 
Laplace- Fourier transforms (Subsection 4.2). Next, we write out the formulas correspond- 
ing to the particular cases N £ {2,3,4} (Subsection 4.3). Finally, we compute the distri- 
bution functions of the vectors (X(t),m(t)) and (X(t),M(t)) (Subsection 4.4) and write 
out the formulas associated with N G {2,3,4} (Subsection 4.5). Although A" is assumed 
to be even, all the formulas obtained in this part when replacing A" by 3 lead to some 
well-known formulas in the literature. 



4.1 Laplace- Fourier transforms 



Theorem 4.1 The Laplace- Fourier transform of the vectors (X(t), M(t)) and (X(t),m(t)) 
are given, for 5ft(A) > 0, \x G R, v > 0, by 



E, 



E, 



f 

J 

I 

Jo 



-Xt+ifxX(t)-vM(t) ^ 



e -\t+inX(t)+vm(t) ^ t 



mv\-(iii-u)9j) n (V\-ifie k y 

j&J keK 



,{i^+v)x 



(4.1) 



n(VA-i^) \[{V\-{iii + v)e k )' 

jeJ k&K 
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Proof. We divide the proof of Theorem 4.1 into four parts. 
• Step 1 

Write functional (3.1) with H x (t) = X(t), K x (t) = M(t) or K x (t) = -m(t) and 
Ix(t) = 1: 

r+oo r+oo 
F+{\,ti,v)= e -\t+i»X(t)-vM(t) dt and F -(X,p, V ) = / e -M+inX(t)+vm(t) dt 

Jo Jo 
So, putting X nyk = X(t Uik ), M n (t) = max ^ t I I1 (s) = max ^^ L2 ™tJ X n>j where [-J 
denotes the floor function, and next M nk = M n (t nk ) = max ^j<fe X n j, (3.2) yields, e.g., 

F+ n (A,/X,I/) = ^ e -A*n,*+i^ n ,fc-I'Af„,fc - 

A;=0 

The functional [\,^,v) is a function of discrete observations of X. Our aim is to 
compute its expectation, that is to compute the expectation of the above series and next 
to take the limit as n goes to infinity. For this, we observe that, using the Markov property, 

k 



E 



-Xtn 



j=0 



< (e 



-$t(\)/2 n \k 



k 

/•••/ e^- ra ^(l/2 n ;x-x 1 )---p(l/2";x fc _ 1 -x fc )dx 1 ---^. 

^ {k + l){pe-^/ 2n ) k . 

So, if 3?(A) > 2 n mp, the series ^ E x [e _A * n -'= +vXn '' i: ~ jyMn ' fc ] is absolutely convergent and 
then we can write the expectation of (X, fi, u): 

1 - e" A /2" 



E x [F+ n (A, M ,f)] 



A 



k=0 



<p+ k (n,v;x) for »(A) > 2™ In/? (4.2) 



with 



-(v-in)M n>k -in(M n>k -X n , k ) 



However, because of the domain of validity of (4.2), we cannot take directly the limit as n 
tends to infinity. Actually, we shall see that this difficulty can be circumvented by using 
sharp results on Dirichlet series. 

• Step 2 

Putting z = e~ A / 2 ™ and noticing that e _A * n - fe = z k , (4.2) writes 

1 -z 



E x [F+ n (X,n,v)] = — ^-^^(/z^x)**. 



k=0 



The generating function appearing in the last displayed equality can be evaluated thanks 
to an extension of Spitzer identity, which we recall below. 

Lemma 4.2 Let ($, k )k>i be a sequence of "i.i.d. random variables" and set Xq = 0, 
X k = Y^=i Cj f or k ^ 1, and M k = maxo^j^fcATj for k ^ 0. The following relationship 
holds for \z\ < 1: 



^E[e^ x "- vMk ] z fe = exp 



k=0 



co k 

sr E ^x k -„x^ - 
k=i 



ii 



Observing that 1 — z = exp[log(l — z)\ = exp[— J2T=i zk Lemma 4.2 yields, for 



E x [F+( A, /*,!/)] 



e (v*-0* exp 

A 



4e 



(4.3) 



where 



V> + (/W) = E 
= E 



^x{t)-vx{ty 



- 1 



e ifiX(t) _ A i 



HX(t)<0} 



+ E, 







(i/x-i/)X(t) 



l) l{X(t)>0} 



= | ^_i) p ( i; _ e)de + jf <X> ^(v-^_i)p( t; _^) de . (4.4) 

We plainly have |Y> + (/x, v;t)\ < 2/?, and then the series in (4.3) defines an analytical function 
on the half-plane {A G C : 5ft(A) > 0}. It is the analytical continuation of the function 
A i — ► E x [i ? ^(A, fi, v)\ which was a priori defined on the half-plane {A G C : 3?(A) > 
2 ra lnp}. As a byproduct, we shall use the same notation E x [.Fj (A, /x, v)\ for 9ft(A) > 0. 
We emphasize that the rhs of (4.3) involves only one observation of the pseudo-processus 
X (while the lhs involves several discrete observations). This important feature of Spitzer 
identity entails the convergence of the series lying in (4.2) with a lighter constraint on the 
domain of validity for A. 

• Step 3 

In order to prove that the functional F% (A, n, u) is admissible, we show that the series 



is absolutely convergent for 5ft(A) > 0. For this, we invoke a 



lemma of Bohr concerning Dirichlet series ([5]). Let aj t e~ /3kX be a Dirichlet series of the 
complex variable A, where (afc)fceN is a sequence of complex numbers and (/3fc)fceN is an 
increasing sequence of positive numbers tending to infinity. Let us denote a c its abscissa 
of convergence, a a its abscissa of absolute convergence and o~b the abscissa of boundedness 
of the analytical continuation of its sum. If the condition limsup^^ \n(k)/0k = is 
fulfilled, then a c = a a = u\,. 

In our situation, we will show that the function of the variable A lying in the rhs in (4.3) 
is bounded on each half-plane K(A) ^ e for any e > 0. We write it as 



exp 



-At 



n,k 



k=i 



exp 



00 —\t 



E 



■E, 



jliX nk 



x exp 



E 

k=i 



Xtjr, 



■E, 



(e^~^ - l) l {Xn ^o } 



For any a G C such that 5ft(a) ^ 0, we have 



Ei 



(e"* W - l) l {xm0} 



E 

+oo 



; «tV"x(i) _ -A j 







1{X(1)^0} 



x|p(l;-0|de 



< 2f ? |a|t 1/Af 



where we set £ = / + °° £ |p(l; — 01 ^£ < +oo) and we used the elementary inequality 
|e^ — 1| ^ 2|£| which holds for any (6 C such that 9i(C) ^ 0. Similarly, 



E ( 



^x(t) 



l) !{x(t)<o} 



< 2 e |a|t 1 /". 
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Therefore, 



k=i 



k 



JsO (or < 0)} 



< 2£?|«| ^ 



71 A_ , I 



2- — + 

k=l k=l L n,k 



J„ 1 " ta 



k n ' k 

/-tn.fc+i e -SR(A)t 



2g\a\ f 
Jo 



l/N 
+00 e -SR(A)t 



fl-l/JV 



2r(i/JV)g|a| 



This proves that the rhs of (4.3) is bounded on each half-plane R(A) ^ e for any e > 0. 
So, the convergence of the series lying in (4.2) holds in the domain 3J(A) > and the 
functional F^(X,fj,,v) is admissible. 

• Step 4 

Now, we can pass to the limit when n — > +00 in (4.3) and we obtain 

for »(A) > 0. (4.6) 



E :E [F+(A,^,^)] = ie^-^exp 



A similar formula holds for 



From (4.4), we see that we need to evaluate integrals of the form 

r+00 jj. r+00 

/ e -A* ^ / ( e «£ _ i) p(t; _£) de for R( a ) < 
Jo 1 Jo 

and 

f + °° tit f° 

/ e~ xt — / (e af - l)p(t; -f ) d£ for »(a) > 0. 
We have, for ft(a) < 0, 

/•+00 Jj. >>+oo 

Jo 1 Jo 

r+00 r+00 r+00 
= dt e- ts ds ( e «_i)p(t ; -£)d£ 

J0 JA JO 

r+00 r+00 r+00 

= ds (e^-l)de/ e" 4 * 

JA JO JO 

/• + OO /-+00 / \ 

= / da {e a t -1) £ i d ? P uttin S CT = V^) 

j^ jo v eJ y 



E 



+00 



VA 



— a9j 



(4.7) 



In the last step, we used the fact that the set {0j,j £ J} is invariant by conjugating. 
In the same way, for 3ft(a) ^ 0, 



r+00 j± i-O 

/ e - A *f/ (e*-iMt;-*K = £ 

JO J J-00 



(4.8) 
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Consequently, by choosing a = ifi in (4.7) and a = ifi — v in (4.8), and using (2.1), it 
comes from (4.4): 



cxp 



-oo ^£ 
e~ xt ij+(p,v;t)j 



jeJ keK 



From this and (4.6), we derive the Laplace-Fourier transform of the vector (X(t), M(t)). 
In a similar manner, we can obtain that of (X(t),m(t)). The proof of Theorem 4.1 is now 
completed. ■ 

Remark 4.3 Any of both formulas (4.1) can be deduced from the other one by using a 
symmetry argument. 



dist 



For even integers N, the obvious symmetry property X = — X holds and entails 



E, 



Jfj,X(t)+u mino^s^t X(s) 



= E 
= E 



-i^X(t)-v max ^i X(s) 



Observing that in this case {6k, k G K} = {—6j,j G J}, we have 



and 



TT re 

n re 



n , w 

VX + i^6 k 

n ^ 



which confirms the simple relationship between both expectations (4.1). 

If iV is odd, although this case is not recovered by (4.1), it is interesting to note the 
asymmetry property X + d = —X~ and X~ d = —X + where X + and X~ are the 
pseudo-processes respectively associated with kn = +1 and kn = — 1. This would 
give 



Ei 



JfJ,X+ (t)+u min <g s ^ t X+ (s) 



= E 
= E 



-i/iX (t)+z/ min ^ 3S ; t (-X (a)) 



-j|itX (t)-l'mH ^<; i I (s) 



Observing that now, with similar notations, G J + } = {— k ,k G K } and 

{0£, k G = {— 9~, j G J - }, the following relations hold: 



n 

jeJ+ 



Vx 



and 



n 



n 

fcex- 



n 



Vx 



Vx + i/j.8: 



Vx 



Hence (X + (f), m + (t)) and (X (t),—M (t)) should have identical distributions, 
which would explain the relationship between both expectations (4.1) in this case. 
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Remark 4.4 By choosing v = in (4.1), we obtain the Fourier transform of the A-potential 
of the kernel p. In fact, remarking that 



JV-l 



HiVX-ifiOj) Y[(V\-ifi9 k )= H(V\-ifi9i) = \-K N (iti) N , 

j&J keK 1=0 



(4.1) yields 



E , 



/ 

J o 



e -\t+i»x(t) dt 



tflX 



A — KN(ifl) N 



which can be directly checked according as 



/" 

Jo 



xt E x 


' e ifiX{t) 


dt = 






Jo 



+ 0O 



t ifj,x-(\-K N (in) N )t ^ 



4.2 Density functions 

We are able to invert the Laplace-Fourier transforms (4.1) with respect to \i and v. 



4.2.1 Inverting with respect to v 
Proposition 4.5 We have, for z ^ x, 

P + OO 



/" 

JO 



e~ xt dtE x 



and, for z ^ x, 

f+OO 



/ 

Jo 



e~ M dtE x 



e^ X{t \M{t) G dz jdz = 



e^ X{t \m{t) G dz /dz 



X (l-#J)/N e itMX 

n (Vx-iMt 



Em* 



(ifi-8j 1 Vx)(z-x) 



X (l-#K)/N e i^x 



3&J 



U(V\-in9 

J6J 



fee 



(4.9) 

(•i/j. — 0fe i V / A )(z— 



Proof. Observing that {%,J € J} = {Oj,j € J} = {l/9j,j € J}, we have 
1 1 \-# J l N 



jeJ ikj V 



Applying expansion (2.11) to x = (ip - u)/ V\ yields: 



1 



x -#j/n y- = a (i-#j)/jv y- _ 



(4.10) 



Writing now 

we find that 
o 



f — ifi + 9j a/A A 



+oo 



e (in-6j t V\){z-x) ^ 



jnX{i)-vM(i) 



dt 
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^(1-#J)/N e i^x 

n (Vx-iM 

We can therefore invert the foregoing Laplace transform with respect to v and we get the 
formula (4.9) corresponding the case of the maximum functional. That corresponding in 
the case of the minimum functional is obtained is a similar way. ■ 

Formulas (4.9) will be used further when determining the distributions of (r+ , X(t+)) 
and (t-,X(t-)). 




4.2.2 Inverting with respect to [i 

Theorem 4.6 The Laplace transforms with respect to time t of the joint density of X(t) 
and, respectively, M{t) and m(t), are given, for z ^ x V y, by 

1 



»+oo 

e~ xt dt¥ x {X(t) G dy,M(t) G dz}/dydz = ^ip x (x - z)i/>\(z - y), 



f 

Jo 

and, for z < x Ay, (4-11) 

/ e- xt dtP x {X(t) G dy,m(t) G dz}/dydz = - t(; x (x - z) ip x (z - y), 
Jo A 



r+oo 
10 

where the functions ip x and ip x are defined by (2.6) 



Proof. Let us write the following equality, as in the previous subsubsection (see (4.10)): 



1 



n (Vx-ifid k ) 



_ x (i-#k)/n _ 
keK W 



GkBk 



e k Vx' 



Set 



+oo 



f 

G(X,fi;x,z) = 

Jo 

We get, by (4.9) and (2.1), for z ^ x, 

X (i-#J)/N e n 



e~ xt dtE, 



e^ X{t) ,M(t) G dz /dz 



G(X,n;x,z) 



VMie (iMiv/I)M 



keK 



_ x (2-#J-#K)/N e i fJ ,x 

keK ¥~0kVX 



jeJ,keK 



ifi - 9k VX 



Writing now 

gives 

G(X,n;x,z) 



ifi - 6 k VX 



f 

j — c 



dy 



-X 2 ' N ~ 



J — c 



OjAj 6 k B k e N ^(^-e k y+(e k ~e,)z) 

jeJ,keK 



dy 
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and then 

r+oo 

/ e~ xt dt P x {X(t) G dy, M(t) G dz) jdy dz 
Jo 

j&J,keK 

This proves (4.11) in the case of the maximum functional and the formula corresponding 
to the minimum functional can be proved in a same manner. ■ 

Remark 4.7 Formulas (4.11) contain in particular the Laplace transforms of X(t), M{t) 
and m(t) separately. As a verification, we integrate (4.11) with respect to y and z sepa- 
rately. 

• By integrating with respect to y on [z, +oo) for z ^ x, we get 

r+oo 

/ e~ Xt dt¥ x {m(t) G dz}/dz 
Jo 

= _ A 2/^v-i djAj e~ e ^y-^ dy k B k e **^(*-*) 

jeJ Jz keK 
= _ A ViV-i £ A, k B k e 9 ^-^ 

jeJ keK 

= -A^X;^^ = Ux(x-z). 
keK 

We used the relation J2jeJ A? = see Subsection 2.3. We retrieve the Laplace 
transform (2.5) of the distribution of m(t). 

• Suppose for instance that x ^ y. Let us integrate (4.11) now with respect to z on 
(— oo,x]. This gives 

r+oo 



r+oo 

/ e~ xt dt¥ x {X{t) edy}/dy 
Jo 



jeJMK 



\l/N-l \ ^ fyA? 6>fc-Bfc Jj'Vxix-y) 

^ e k -0j 



i 





jeJ \fce_(Y K J / 
= l A i/jv-ij-^ e ^( a! - w ) ) 

where we used (2.10) in the last step. We retrieve the A-potential (2.3) of the 
pseudo-process {X{t))t^o since 

r+oo r+oo 

/ e~ M dt¥ x {X(t) G dy}/dy = e~ xt p(t;x - y) dt. 

Jo Jo 
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Remark 4.8 Consider the reflected process at its maximum (M(t) — X(t))t^o- The joint 
distribution of (M(t), M{t)-X{t)) writes in terms of the joint distribution of (X(t), M(t)), 
for x = (set P = Po f° r short) and z, ( ^ 0, as: 

F{M(t) G dz, M(t) - X(t) G d(} = F{X(t) G z - d(, M(t) G cfe}. 

Formula (4.11) writes 

r+oo 

/ \e~ xt dtF{M{t) G tfe,M(t) - X(t) G dC}/ctedC = PaWa(-C) 
jo 

r+oo r+oo 

= / A e _A< P{M(t) G dz}/dz x / Ae" Ai dtP{-m(t) G d(}/d(. (4.12) 
jo jo 

By introducing an exponentially distributed time T\ with parameter A which is indepen- 
dent of (X(t)) t ^o, (4.12) reads 

P{M(T A ) G dz, M(T\) - X{T X ) G dQ = P{M(T A ) G dz}F{-m(T x ) G d(}. 

This relationship may be interpreted by saying that —m(T\) and M(T\) — X(T\) admit 
the same distribution and M(T\) and M(T\) — X(T\) are independent. ■ 



Remark 4.9 The similarity between both formulas (4.11) may be explained by invoking 
a "duality" argument. In effect, the dual pseudo-process (X*(t))t^o of (X(t))t^o is defined 
by X*(t) = —X(t) for all t ^ and we have the following equality related to the inversion 
of the extremities (see [12]): 

F x {X(t) G dy, M(t) G dz} jdy dz = F y {X* (t) G dx, m* (t) G dz} jdx dz 

= F- y {X(t) G d{-x),m(t) G d(-z)}/dxdz. 



Remark 4.10 Let us expand the function ip\ as A — > + : 



^(o = ^Em4e 1i ^t^ + °( a(#j " 1)/jv ) 

ieJ L z=o 



#J_1 / \ \(l+l)/N(l 

z=o \ jeJ / 

We have by (2.11) (for a: = 0) £ jGi/ 0j +1 A,- = for < I < # J - 2 and T, je jdf J Aj = 
(-1) #J-1 IW^ Hence 



(#«/-!)! 



Similarly 



mo - + (-i)* k n ^ -* K ~ 1 - x#k/n - 



(4.13) 



(4.14) 



As a result, putting (4.13) and (4.14) into (4.11) and using (2.1) and UiLq 1 °i = (-1)^" W 
lead to 



r+oo 

/ e~ xt dtF x {X(t) G dy,M(t) G dz}/dydz ~ kjv 

Jo A^0+ 



( x _ z )#J-l( z _ y )#^-l 

a-o+ " W (#J-l)!(#tf-l)! ' 
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By integrating this asymptotic with respect to z, we derive the value of the so-called 
O-potential of the absorbed pseudo-process (see [19] for the definition of several kinds of 
absorbed or killed pseudo-processes): 

r+oo ra i _ \#J-lt _ \#K-1 

I P. Wt ) 6 dy ,M(t) 4 a }/d y = HPV»/ LL_1_J[_ dz . 



4.2.3 Inverting with respect to A 

Formulas (4.11) may be inverted with respect to A and an expression by means of the 
successive derivatives of the kernel p may be obtained for the densities of (X(t),M(t)) 
and (X(t),m(t)). However, the computations and the results are cumbersome and we 
prefer to perform them in the case of the distribution functions. They are exhibited in 
Subsection 4.4. 

4.3 Density functions: particular cases 

In this subsection, we pay attention to the cases N G {2,3,4}. Although our results are 
not justified when N is odd, we nevertheless retrieve well-known results in the literature 
related to the case N = 3. In order to lighten the notations, we set for, 9J(A) > 0, 

r+oo 

&\(x,y,z) = / e~ xt dtF x {X(t) €dy,M(t)€dz}/dydz, 
Jo 

r+oo 

V\(x,y,z) = / e~ xt dtF x {X(t) G dy,m(t) G dz}/dydz. 
Jo 

Example 4.11 Case N = 2: using the numerical results of Example 2.1 gives 
<Px(0 = y/\e VX * and ^ x (0 = VXe"^, 

and then 

*x(x,y,z) = e^+v-V l {z > xVy} and * x (x,y,z) = e ^ 2 — ») t {z ^ xAy} . 



Example 4.12 Case N = 3: referring to Example 2.2, we have 

• for K3 = +1: 

MO = ^Ae^, 

MO = -^(e "*-e ^) = ^ r e >*8m{—V\t), 

which gives 

<Mw) = -^ e ^+^-|^) s in^^A(z-y))l { ^ Vy} , 
*x(x,y,z) = -^e^h-k*-v) s in(^n(x-z^l { ^ xAy y, 
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MO = 

MO = V\e~^, 
which gives 



< i sin | ^ \ . 



2 



z - x) 1 



Example 4.13 Case N = 4: the numerical results of Example 2.3 yield 



MO 
MO 

which gives 
®\{x,y,z) 



V2 



V2 



sin I § 



1 e ^(x+y-2*) 



Va 
1 

7x 



1 



cos ^-^= (x - y)j - cos (x + y- 2z) 

cos I — ^= {x-y)\ - cos I — (x + y - 2z) 



4.4 Distribution functions 

In this part, we integrate (4.11) in view to get the distribution function of the vector 
(X(t),M(t)): P x {X(t) < y,M(t) < z}. Obviously, if x > z, this quantity vanishes. 
Suppose now x ^ z. We must consider the cases y ^ z and z ^ y. In the latter, we 
have F x {X(t) ^ y,M(i) ^ z) = P{M(i) < z} and this quantity is given by (2.7). So, we 
assume that z ^ x V y. Actually, the quantity ¥ x {X(t) ^ y, M(t) ^ z} is easier to derive. 

4.4.1 Laplace transform 



Put for 31(A) > 0: 



h(x,y,z) 



r+oc 

/ e- xt ¥ x {X(t)^y,M(t)^z}dt 
Jo 

r+oc 

/ e- xt ¥ x {X(t)^y,M(t)^z}dt. 
Jo 



The functions F\ and F\ are related together through 

F x (x, y, z) + F x (x, y, z) = f(A; x - y) 
20 



(4.15) 



where ^ is given by (2.4). Using (4.11), we get 

F\(x,y,z) 

/y r-\-co r-\-oo 
/ / e- M P x {X{t)€d£,M(t)€dC}dt 
-co J z JO 

We plainly have ( ^ z ^ xV y ^ xV £ over the integration set (—00, y] x [z, +00). So, the 
indicator l{£^ x v£} i s useless and we obtain the following expression for F\. 



Proposition 4.14 We have for z ^ xV y and 9i(A) > 0: 

/•+00 1 
J e- xt F x {X(t)^y^z^M(t)}dt = - £ 



j&J,k£K 



9j - 9k 



and for z ^ x Ay: 

[ +0 °e- xt F x {X(t)^ y >z^m(t)}dt = ± V ^^--H^G*-*). 

As a result, combining the above formulas with (4.15), the distribution function of the 
couple (-X"(i),M(i)) emerges and that of (X(t),m(t)) is obtained in a similar way. 

Theorem 4.15 The distribution functions of (X (t) , M (t)) and (X(t),m(t)) are respec- 
tively determined through their Laplace transforms with respect to t by 
r+00 



r+co 

/ e- xt F x {X(t) <j/,M(t) ^ z} dt 
Jo 



= < 



( 1 x - OjAjBk ej V\(x-z)+e k V\{z- y ) J_ y- e h V\{x- 
X 2^ e Q e iVA ^ 

j&J,k&K J k&K 



j&J,k£K 



y) ify^x^z, 
(4.16) 
ifx^y^z, 



and 



r+00 

/ e- xt P x {X(t)^y,m(t)^z}dt 
Jo 

'1 ^ A^ fc g fc e ^ (3! -,) +0fc ^-v) + J- ^ e ^(»-y) i/_<*<y, 

/\ _____ C7 •? 1/ fa - V /\ 



= < 



j£j,keK 



j&J,k£K 



if z ^y ^x. 



4.4.2 Inverting the Laplace transform 

Theorem 4.16 T/ie distribution function of (X (t) , M (t)) admits the following represen- 
tation: 

F x {X(t)^y^z^M(t)} = £ a fcm f T ^(a- x- z) ds da (4.17) 
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where Iko is given by (5.14) and 

NB k AjOLj r 



the ctj m 's being some coefficients given by (4-18). 

Proof. We intend to invert the Laplace transform (4.16). For this, we interpret both 
exponentials e 9j ^( x ~ z ) an d e® k ^( z ~v) as Laplace transforms in two different manners: 
one is the Laplace transform of a combination of the successive derivatives of the kernel p, 
the other one is the Laplace transform of a function which is closely related to the density 
of some stable distribution. More explicitly, we proceed as follows. 

• On one hand, we start from the A-potential (2.3) that we shall call <E»: 

jeJ 

Differentiating this potential (# J — 1) times with respect to £ leads to the Vander- 
monde system of # J equations where the exponentials e 6 ^ are unknown: 

£0*+V^ = NX i-H+i)/N |^(A;0 for < I < # J - 1. 

jeJ 

Introducing the solutions ctj m of the # J elementary Vandermonde systems (indexed 
by m varying from to # J — 1): 

we extract 

e -ijw = n*t a ^ g^(A-o 

A ^ \(m+l)/N dx mK ^ } 

m=0 



f + °° -xt, f f' # y Na jm d m p ds 

Jo e Jo £f r(^i) <9* m ls ' u (t - s )i-(-+D/^v • 

The explicit expression of ctj m is 

lli 6 j\{j}W -^j) llzeJ^ 

where the coefficients Cj q , ^ q ^ # J — 1, are the elementary symmetric functions 
of the 0;'s, I G J \ {j}, that is Cjo = 1 and for 1 ^ q ^ #J — 1, 

c j9 = a,(^,ZG J\{j})= Yl <V"<V 

h,...,i q eJ\{j} 

l 1 <--<lq 

On the other hand, using the Bromwich formula, the function £ i — > e Sk can be 
written as a Laplace transform. Indeed, referring to Section 5.2.2, we have for k £ K 
and £ ^ 0, 

n/- r +°° 
J k V\( = I e ~ xt I k0 (t;t)dt 



10 

where I^q is given by (5.14). 
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Consequently, the sum lying in Proposition 4.14 may be written as a Laplace transform 
which gives the representation (4.17) for the the distribution function of (X(t),M(t)). ■ 



Remark 4.17 A similar expression obtained by exchanging the roles of the indices j and 
k in the above discussion and slightly changing the coefficient a^ m into another bj n may 
be derived: 

F x{ x(t)^y^z^M(t)}= E ^ f Q [ (*; * - v) l°[ S ~/-(n^)/N ds d ° ( 4 - 19 ) 



where 



°0n ~ r(n +l\ Z_> 



However, the foregoing result involves the same number of integrals as that displayed in 
Theorem 4.16. ■ 



4.5 Distribution functions: particular cases 

Here, we write out (4.16) and (4.17) or (4.19) in the cases N E {2,3,4} with the same 
remark about the case N = 3 already mentioned at the beginning of Subsection 4.3. The 
expressions are rather simple and remarkable. 

Example 4.18 Case N = 2: the double sum lying in (4.16) reads 



E 

j£j,k£K 



9k — Oj 



with = -i, and then 



F x (x,y,z) = < 



1 

2A 

1 1 

A ~ 2A L 



e V\(x-y) + e ^/X(x+y~2z) 



if y < x z, 
if x ^ y ^ z. 



Formula (4.17) writes 

rt PS 

F x {X(t) ^ y < z < M(t)} = a 00 J J p(<r\ x - z) I Q0 (s -<r;z-y) -j= 



with 



p(t;0 



l 



TTt 

The reciprocal relations, which are valid for £ ^ 0, 



e 4t . 



2^/\ 



and = 2v / Aai ^(A;O 



imply that ctio = 1. Then aoo = r(i/2) = 7/?' ^ n ^ e °th er hand, we have for £ ^ 0, 

by (5.14), 



/oo(*;0 = 



2itt 



■/ 

t r /•+oo /"+oo 
27rt [J J 

2vrt J_ 



dX 



20Ft 3 / 2 



e 4t . 
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Consequently, 



2 

Using the substitution a = together with a well-known integral related to the modified 
Bessel function K 1 / 2 , we get 



(x-z) 2 (z-y) 2 



l' S e 4<T 4(s-ct) g 4s 



00 (^-f) 2 u ^-y) 2 du 

g J2 Au 



u; 



3/2 



^/VT (2z-x-y) 2 

e 4 3 



Then 



,{X(t) ^y^z^ M(t)} = — 



(2z-x-yy 

ire 4s 



ds. 



27r Jo y/a(t - s) 

Finally, it can be easily checked, by using the Laplace transform, that 



. _ (2z-x—y) 
" " 4s 



/ - ds = ^ dt = 2ir p(t; -0 dt. 

Jo ycryt — s) J2 Z -x~ y yt J2 Z -x- y 

As a result, we retrieve the famous reflection principle for Brownian motion: 
F x {X(t) ^ y s$ z M(t)} = F{X(t) ^2z-x-y}. 



Example 4.19 Case N = 3: we have to cases to distinguish. 



Case K3 = +1: the sum of interest in (4.16) reads here 

02 — &0 



with 



Bi 

9i-e 



h-0o 
~- sBk = ~h and then 



F\{x,y,z) = < 



2 

3A 



v. yS> cos (-^- \^(x — y) 



if y ^ x ^ z, 



1 1 

A ~ 3A 



, \f\{x-y) 



V3 



z y 2 z ) cos ( — V\(z - y) 



if x ^ y ^ z. 



We retrieve the results (2.2) of [3]. Now, formula (4.17) writes 

F x {X(t) ^y < z < M(t)} = I I p{a;x-z)(a w I w +a 2 oho){s-a-z-y) 

Jo Jo 

where, by (2.2), 

i r+oo 

p(t;0 = - cos(£A - U 3 ) dX. 
n Jo 



ds da 

(t - s) 2 /3 
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The reciprocal relations, for £ ^ 0, 



imply that aoo = 1- Then 

aw = 

O20 = 

Consequently, 

F x {X(t) ^ y < z < M(t)} = 
with, for £ ^ 0, by (5.14), 

g(t;0 = (/io + / 20 )(i;0 



3Si A)aoo 



1 



r(i/3)e -(9i r(i/3)' 

3B 2 A a 00 1_ 

r(i/3) e -e 2 ~ r(i/3)' 



p(<r; x - z) g(s - cr; z - y) 



r(i/3) 



JO 



(t - s)2/3 



2vrt 



Jo Jo 

+oo 



27Tt 
JO 



Jo 



r+oo 



+ 00 i7r . f+OO in 



7T 



£A + - dA. 



Case K3 = —1: the sum of interest in (4.16) reads here 

01 — 0Q 01 — 02 



with 



9o A Q 

01-00 



-± e i7r / 3 and |^ = -± e" i7r / 3 , and then 



F\{x,y,z) 



1 

3A 

if y < x ^ z, 



7T 



1 1 

A ~ 3A 



2e^-^cos V\(x-y) 



+ e 2 (2 x+ y 2 z ) cos 
if x ^ y ^ z. 

We retrieve the results (2.2) of [3]. Next, formula (4.19) writes 
¥ x {X(t) < y ^ z M(t)} 

X] 6 i0 / / P(cr;z-y)l j0 (s-a;x- z 

Wml JO JO 



(is tier 



ie{0,2} 

ft /-S 



= / / P(v,z-y) (booho + ^20^20) (s - a; x - 
Jo Jo 



(t - a)2/3 

(is dcr 



(t - s)2/3 
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where, by (2.2), 



I r+oc 

pit- = - cos(£A + tX 3 ) dX. 
* Jo 



From the reciprocal relations, which are valid for £ ^ 0, 

$ ( A ^) = ^7T and -e-^ = 3A 2 / 3 /3 10 <I>(A;e) 
we extract the value (3io = — 1. Therefore, 



2l7T 



&00 
b 2 

Consequently, 



30 O A O 0iBi/?io e— 



r(i/3) 0i-0 o r(i/3) 



■>ITT 



30 2 A 2 exBrfw es 



r(i/3) 0i-0 2 r(i/3) 



1 /* ^ d d 



r(i/3) JJ ^ ' y ' yv ' '(t- s )2/3 



where, for £ ^ 0, by (5.14), 



r>+00 o j,r ._ /"+00 

-0n / 

2vrf 



L Jo Jo 

+ 2 e"f / e -^+e 2 eir 5 A dX + d J e -tX^e~^X dX 
Jo Jo 

V3 £°° e-' A3+ « A dX + ^ + °° e - tA3 -5 ^ S in £A + |) dA 



1 



Example 4.20 Case N = 4: in this case, we have 
jeJ,keK K 3 

6q — 2 0\— #2 

I 9 3 A 3 B e e 3 &\(x-z)+Oo V\(z-y) + 6 3 A 3 B 1 & 8 3 ^X(x-z)+9 1 tf\(z-y) 
6q — 3 0i — 03 

with 

2 A 2 B _ i 3 A 3 B 1 _i_ 2 A 2 B 1 _ e~^/ 4 3 A 3 B o _ e^/ 4 
o -0 2 4' 0i- 03 4' 0i- 2 2>/2 ' O o -0 3 2V2' 
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Hence, 



F\{x,y,z) 



1 

2A 



-$g(x-y) ( ^ i \ M(x+y-2z) 

^ y y) cos I — y=(x - y) J - e y ' 



V2 

if y < x ^ z, 



x cos —i={x — y) J — sin — ^(a? — y) J — sin —^=(x + y — 2z) 



V2 



VV2 



1 1 

A ~ 2A 



is 

e V5 



COS 



^ (I _ v) ) +e -£<«»-*> 



x I cos I —i={x — y) J — sin — ^(a? — y) J — sin — ^(a? + y — 2z) 



^2 ) 
k if x ^ y ^ 2. 

We retrieve the results (3.2) of [3]. Now, formula (4.17) writes 

ft PS 



\V2 



F x {X(t) ^ y ^ z < M(i)} = / / p(o-;x - z) (aoo-foo + aio^io)(s - o-;^ - y) 

JO JO 

+ yy — (cr;x-z)(a i/oo + anAo)(s-o-;z-y) 
where, by (2.2), 

i r +co 

p(t;0 = - e~ tx co S (£\)d\. 

n Jo 

Let us consider the system 

' e 2 e e ^ + 6 z e e ^ = 4A 3 / 4 $(A;0 



ds da 



(t - s) 3 / 4 

ds da 



which can be conversely written 



63 — 02 
4 



# 3 A 3 / 4 *(A;0- >/A^(A;0 
-^ 2 A 3 / 4 cI>(A;0 + Va||(A;0) 



#3 — #2 

or, by means of the coefficients 020, «2i, 030, «3i, 

' = 4(a 2 oA 3/4 «J>(A;0+a2iVAH(A;e)), 

3e e 3 m = 4 (a 30 A 3 / 4 $(A;O + « 3 i>/a||(A;o). 
Identifying the two above systems yields the coefficients we are looking for: 

#3 



«20 

"30 
«21 
"31 



#3 — #2 
#2 



A 2 = 

= A 3 



e * 



>/2 ' 

17V 

e * 

= 7I' 



# 3 — 62 
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and next: 



aw 



an 



JBo_ 

r(i/4) 

4Bi 

rW4) 
r(i/2) 

4Bi 
f(T72) 



^ 2 «20 -43"30 



#2 — 6*0 #3 ~ #0 
^ 2 «20 -43"30 



#2 — 6*1 #3 ~ 0\ 

A 2 a 2 i A 3 a 31 



02 — 6>o #3 — #0 
^2«21 -43031 



v / 2r(i/4)' 
1 

v / 2r(i/4)' 

e 4 



92 — 6*i #3 — #1 



'27T' 

e 4 



2vr 



Consequently, 



; {I(^i/<^M(t)} =11 p(v,x-z)q 1 (s-a;z-y) 

Jo Jo 



ds da 



ft dp 

I I — (a;x-z)q 2 (s-a;z-y) 



(t - S) 3 / 4 

ds da 



Jo 



dx 



with, for £ > 0, by (5.14), 



Qi(t;Q = (a 00 Ioo + a w I 1Q )(t;^) 



\/2r(l/4) 



Uoo + /io)(t;6 



and 



92(*;£) = (aoiioo + an/io) (*;£) = -7= (e 4 ^00 + e * iio)(*; O 

V 27T 



= ^ (ioo + Ao)(*;0 - ^ 0*oo - ho)(t;0- 
Let us evaluate the intermediate quantities (loo ± ho)(t',0 : 



{ioo + homo = ^ 



and 



Then 



i e 



r+00 i7r r+00 ijr 

Jo Jo 

-00 



/•+00 

Jo Jo 

£ f+oo 
TTt ./I) 



/ e"' A cos(£A)dA 

JO 



(/oo-/io)(*;0 



f+°° 4 r+00 
/ e -* A4 -^dA- / e-' A sin(eA)dA 
Jo Jo 



t r+00 

t r+00 . . 

92(*;0 = ^2- t J ^ (cos(eA)+sin(eA)- e -« A ) dA. 
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4.6 Boundary value problem 



In this part, we show that the function x i — ► F\(x, y, z) solves a boundary value problem 

d N 
IP 



related to the differential operator V x = kn -ttt- Fix y < z and set F{x) = F\(x, y, z) for 



x G (—00, z\. 

Proposition 4.21 The function F satisfies the differential equation 

I XF(x) — 1 for x £ (—00, y), 

[ XF(x) forx € (y,z), 

together with the conditions 

F {l) ( Z -) =0 /or0<i<#^-l, (4.21) 

i?(0( y +) _i?(0( y -) =0 /orO<KiV-l. (4.22) 

Proof. The differential equation (4.20) is readily obtained by differentiating (4.16) with 
respect to x. Let us derive the boundary condition (4.21): 

F (0( Z -) = y Q T A i B * c e„ N Vx(z-«) , * /N y e i + i c e„ N Vx(z-v) 



x i/n-i ^2 



k&K 



yh A] Bk + d -±- 



e 



where the last equality comes from (2.11) with x = 6 k which yields y = — Jf 

^0 k - 0j NB k 

Condition (4.22) is quite easy to check. 



e\ +1 A, el +1 



Remark 4.22 Condition (4.22) says that the function F is regular up to the order n — 1. 
It can also be easily seen that F^ N \y + ) — F^ N \y~) = kn which says that the function 
has a jump at point y. On the other hand, the boundary value problem (4.20)- 
(4.21)-(4.22) (the differential equation together with the N + ffj conditions) augmented 
of a boundedness condition on (— 00, y) may be directly solved by using Vandermonde 
determinants. ■ 



5 Distributions of (t+,X(t+)) and [r a ,X(r a )) 

The integer N is again assumed to be even. Recall we set r+ = inf{i ^ : X(t) > a} 
and t~ = inf{t ^ : X(t) < a}. The aim of this section is to derive the distributions 
of the vectors (t+,X(t+)) and (t~,X(t~)). For this, we proceed in three steps: we first 
compute the Laplace-Fourier transform of, e.g., (t+,X(t+)) (Subsection 5.1); we next 
invert the Fourier transform (with respect to fi, Subsubsection 5.2.1) and we finally invert 
the Laplace transform (with respect to A, Subsubsection 5.2.2). We have especially ob- 
tained a remarkable formula for the densities of X{t^) and X(t~) by means of multipoles 
(Subsection 5.4). 
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5.1 Laplace- Fourier transforms 



We have a relationship between the distributions of (t£,X(Tq)) and (X(t),M(t)), and 
between those of (t~ ,X(t~)) and (X(t),m(t)). 



Lemma 5.1 The Laplace- Fourier transforms of the vectors (t£ ,X(t£)) and (r a , X(r a )) 
are related to the distributions of the vectors (X(t),M(t)) and (X(t),m(t)) according as, 
for K(A) > and 



E, 



E, 



D -\T++inX(T+) 



-Xr a +i^X(r a ) 



{x-K N (ifi) N ) r 

Jo 

(\- KN (i») N ) r 

J 



e" At E. 



e~ xt E x 



e^ x{t \M(t) >a dt forx^ a, 



(5.1) 



e^ XW ,m(t) < a dt for x ^ a. 



Proof. We divide the proof of Lemma 5.1 into five steps. 
• Step 1 

For the step-process (X n (t))t^o, the corresponding first hitting time r+ n is the instant 
t n ,k with k such that X(t n j) ^ a for all j £ {0, . . . , k— 1} and X(t n ^) > a, or, equivalently, 
such that M nife _i < a and M„ )fc > a where M„ ife = max ^ fe X n)j and X nifc = X(t njk ) for 
^ and M n _\ = -co. We have, for x ^ a, 



p -Ar a '; n +i^X n (r Q f; n ) -\t nk -\-i[xX nk 11 

e _ 2^ JL {M n , fc _ 1 <a<M„ ifc } 



fc=0 

00 

£< 

fc=0 



l{M„ ifc >a} !{M»i-i>«} 



(5.2) 



Let us apply classical Abel's identity to sum (5.2). This yields, since 1{m„ _i>a} = an d 
lim^+oo e~ M ^+^ x n, k l {Mn k>a} = 0, for R(A) > 0: 



fc=0 



{M n>fc >a}- 



The functional e Ar ^™+ l ^n(f a ';n) i s a function of discrete observations of X. 
• Step 2 

In order to evaluate the expectation of the foregoing functional, we need to check that 
the series 



y * E a; ^e~^*"' fc " t "^^ n ' fc — e~ xtn ' k+1+ifJ/Xn ' k+1 ^j 1 



1 {M njfc >a} 



fc=0 



is absolutely convergent. For this, we use the Markov property and derive the following 
estimate: 



E 



„-^tn.k+ilJ-X n ,fe1| 



c JL {A n ,i^a,...,A nifc ^a} 



< le- A / 2 T 



/a />a 
OO J — DO 



e^ Xk p{l/2 n ;x - Xl ) ■ ■ -p(l/2 n ; x fc _! - x fc ) dxi • • • dx fc 



< (pe-^)/ 2 ")^ 
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We recall that in the last inequality p = f^\p(t;z)\dz < +00. Similar computa- 
tions yield the inequality |E x [e _A * n ''= +JAtXn ' fe ] | ^ (pe~^( x ^ 2 ") k . Because of the identity 
l {M n , k >a} = 1 - l{Afn, fc ^}, we plainly get 



E. 



Upon adding one integral more in the above discussion, it is easily seen that 



As a result, when choosing A such that 5ft(A) > 2 n lnp, we have 



00 

\ X \ J \"^n,k> a ! 



k=0 

• Step 3 



2(1 + pe 



-5R(A)/2" 



1_ pe -Si(A)/2" 



< +OO. 



Therefore, we can evaluate the expectation of e Ar oU+*M^n(TaU) _ gy M ar k ov prop- 
erty we get, for 5ft(A) > 2" hip, 



E 



OO 



e~ xtn ' k E. 



k=0 

00 



-A/2™ „j/j(X nifc+1 — X Uj 



fc=0 



Since Ex^e^"- 1 '^) = e^M^/ 2 " we obtain, for R(A) > 2 n ln/>, 



E, 



= 2 n (^l - e"( A_Kjv( ^ )JV ) /2n ) 



e 1H, X e Jl{M n , fc >a} 



(5.3) 



fc=0 



• Step 4 



In order to take the limit of (5.3) as n tends to 00, we have to check the validity of (5.3) 
for any A such that 9i(A) > 0. For this, we first consider its Laplace transform with respect 
to a: 



f 

J X 



e~ ua E, 



, - At +„ +iftX„ (t+„ ) 



da 



da. 



k=0 

The sum in the above equality writes, using the settings of Subsection 4.1, 

c+00 



e~ va E, 



ifj,X n k -a 
e 1 {M n , fc >a} 



da 
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C M n k 

k=o L Jx 

= Yj [e~ ux E x (e^ Xn ' k ) - E x (e^ Xn ' k - uMn ' k )] 



k=0 
1 
V 

1 

V 



fc=0 fc=0 
e (ifi-v)x ^ 

j_ _ e -(A- Kjv MiV)/2« ~ l_ e -A/2» E 4 F X n ( A ^^)) 



V 

We then obtain 

»+oo 



p '* L n,k 



1 _ e -{X-K N {iti) N )/V- 



~ 1 _ e -A/ 2 ^ exp I ^ £ -T— W 



fc=i 



Jx 



e~ ua E, 



e 

(ifi-v)x 
V 



A-r+ n +i/iX n (r+„) 



da 



1 _ e -(A-K JV M JV )/2" 
1 i _ e -A/ 2 n eX P 



1 00 — \t 

1 v-^ e 



Z> ~^ ^ u > tn > k > 



271 k=i tn ' k 



Inverting the Laplace transform yields, noting that the function a 1 — ► E x 
is right-continuous, 

r<c+ioo 



-Ar+ n +i/iX„(r ( ^ rl ) 



E, 



-Ar tl , ; n +i/iX„(r £ ^ rl ) 



lim / 

2m e ^o+ y c 



D (i/j,-v)x+u(a+e) 



c— too 

JV-, ; n 



1 _ g-(A-Kjv(i/xr)/2 



1 



3 — Ai r , 



exp ™ £ ~i — ^ + (^' 1/5 



2 ™ ifci tn ' fe 



Putting 
with 



1 - e-V 2n 

1/; i) = i) + ^ 2 (i^ - v\ t) 



dv 
v 



(e*m _ ^ i {x(t)<0} ] , ^( Q;t ) = E [(e aX U _ 1) l {x(t)>0} 



the exponential within the last displayed integral writes 

exp w Z> ~ 1 — ^ + ^> ^ ; tra - fc ) 



2 " fel tn > fe 



1 00 _ \+ , 



= exp ( w £ "1 *">*) exp ( w £ "1 ^ 2< ^ ~ 1/5 

V fc=i ™' fc / V fc=i n,fe 

Noticing that 



1 i-c+ioo j 

2ivr ' 



we get 



-Ar+ n +i^X n (T+„) 



1 - 



1 _ e -(\-K N (ili) Ir )/Z» 



GXP W 2^ "1 ^ W tn ' fc ^ 



1 r c+ 

2«vr y c _j, 



c+ioo 



e {a ~ x)v exp 



1 e~^n,k 



dv 
v 
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By imitating the method used by Nishioka (Appendix in [18]) for deriving subtil ex- 
timates, it may be seen that this last expression is bounded over the half-plane 3?(A) ^ e 
for any e > 0. Hence, as in the proof of the validity of (4.2) for 3J(A) > 0, we see that (5.3) 
is also valid for 9J(A) > 0. It follows that the functional e~ Xl ~ a +i v x ( T a ) i s admissible. 

• Step 5 

Now, we can let n tend to +oo in (5.3). For 5ft(A) > 0, we obviously have 
lim 2 n (l- e -(*-«»M N )/*>) =x _ Kn{¥) n, 

and we finally obtain the relationship (5.1) corresponding to r+. The proof of that corre- 
sponding to t~ is quite similar. ■ 

Theorem 5.2 The Laplace- Fourier transforms of the vectors (t+ , X(t+)) and (r~, X(t~)) 
are determined, for §f{(A) > and fry 



E, 



E. 



-At++^X(t-+) 



-Ar a +i^X(-r a ) 



E- 4 J n i 1 

jeJ ie.J\{j} 



n i 1 

fceK JeK"\{fc} 



0, ) e 8k v^^-") e ^ a /or x ^ a. 



(5.4) 



Proof. Using (4.9) gives 

+oo 

e~ At E x 



/" 

JO 



e^ x W,M(t) > a 



(it 



X (l-#J)/N e i^x 

n i v'A ;,/«,.,. 



/•+00 

E 9 ^ / 



(5.5) 



Plugging the following equality 

N-l 

X- KN (ifi) N = HiVX-i^) = YliVX-iiidj) x HiVx-i^k) 

1=0 j&J keK 

into (5.5) and remarking that the set {0j,j € J} is invariant by conjugating yield 



r+oo 

L 1 



e~ xt E 



, M (t) > a 



dt 



= X (l-#J)/N e i^x 3*J ^ ■ 



A, 



-(iix—6j 1 \/\)(x—a) 



tflX 



wvw jeJ ieJ\{j} 



01 e 



-(ifi~9j 1 \/X)(x—a) 



(5.6) 



Consequently, by putting (5.6) into (5.1), we obtain (5.4). 



Remark 5.3 Choosing \i = in (5.4) supplies the Laplace transforms of r+ and r a : 



E, 



E, 



-Ar a 



keK 
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Remark 5.4 An alternative method for deriving the distribution of (t+, A(r+)) consists 
of computing the joint distribution of (X(t), l(_ oo a )(M(t))) instead of that of (X(t),M(t)) 
and next to invert a certain Fourier transform. This way was employed by Nishioka [18] 
in the case N = 4 and may be applied to the general case mutatis mutandis, m 

Remark 5.5 The following relationship issued from fluctuation theory holds for Levy 
processes: if x ^ o, 



E 



-\T a +iflX(T a ) 



/+~ e -^E x [e i ^(*)-«),M(t) > a] dt 
J +00 e^Eo[e^ M W](it ' 



(5.7) 



Let us check that (5.7) also holds, at least formally, for the pseudo-process X. We have, 
by (2.5), 



+oo 



e~ A *E 



;V (M(t)- fl ) )M ^ ^ a 



dt 



r+oo r+oo 

= / e ifl{z - a) dz e- xt dt¥ x {M(t) £ dz}/dz 

J a Jo 



For x = a, this yields, by (2.11), 

P+OO 



/ 

J 



e~ A< E 



MM 



As a result, by plugging (5.8) and (5.9) into (5.7), we retrieve (5.4). 



(5.8) 



(5.9) 



Example 5.6 Case N = 2: we simply have 



E, 



-Xt++i^X(t+) 



= e i^+V\(x-a) for x ^ flj 



Example 5.7 Case N = 3: 



• In the case K3 = +1, we have, for x ^ a, 

ijj_ 



j&J ieJ\{j} ^ 



and, for x ^ a, 

n (i--^ft)^ (x - o) 
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1 :e -|^A(x-a) 



Therefore, (5.4) writes 



E 



3 -Ar++i At X(T+) 



-Ar a +i^X(r a ) 



2 



i{xa— \ \f\(x-a) 



cos I ^^(x - a) - ^ 



for re > a. 



In the case K3 = —1, we similarly have that 

2 



-Ar++i M X(r+) 



\/A(x— a) 



cos I ^^(x - a) + ^ 



E e _Ara + i / i - x '( r o ) 



for x ^ a, 



Example 5.8 Case iV = 4: we have, for x ^ a, 



e^- n 



and, for x ^ a, 



1 e -^^A(x-a) 



II (l-^^)e^(- ) 
fcex ieic\{i:} ^ ^ 



+ e 
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Therefore, (5.4) becomes 



E, 



-Xr++i^X(T+) 



1 



< '< is | —= \f \ (x — a) + — 

V2 4 



+ -77= sin ( —= \f\ (x — a 



for x ^ a, 



-Xr a +ifiX(r a ) 



7T 



cos f — = V^A (x — a) 
1 ^2 4 



+ fr in (^ (x - a) 



for x > a. 



We retrieve formula (8.3) of [18]. ■ 
5.2 Density functions 

We invert the Laplace- Fourier transform (5.4). For this, we proceed in two stages: we first 
invert the Fourier transform with respect to \x and next invert the Laplace transform with 
respect to A. 

5.2.1 Inverting with respect to /x 

Let us expand the product n«eJ\{j} 0- ~ lX ) as 



#J-i 

n (i - 9ix) = £ c jq (- 

lGJ\{j} 9=0 



(5.10) 



where the coefficients Cj q , ^ q ^ # J — 1, are the elementary symmetric functions of the 
6*/'s, I 6 J \ {j}, that is, more explicitly, Cjo = 1 and for 1 < q ^ # J — 1, 



c jq = a q (6i, I G J\{j}) 



E ^ 



i 1 ,...,! 9 gj\{j} 
Il< — <l, 



In a similar way, we also introduce dko = 1 and for 1 ^ q ^ #.K" — 1, 



d kq = <y q (9 t ,leK\{k}) 



E 



l lr ..,l,€K\{fe} 

i!<---<ig 



By applying expansion (5.10) to x = ifi/\/X, we see that (5.4) can be rewritten as 

= E" 4 -'' E ^ n 
jeJ q=0 

#J-1 



E, 



= -Ar++i At X(r+) 



9 



Now, observe that (—ifi) q e lfia is nothing but the Fourier transform of the q th derivative 
of the Dirac distribution viewed as a tempered Schwartz distribution: 

/+0O 
e ilM Si q \z)dz. (5.11) 
-oo 

Hence, we have obtained the following intermediate result for the distribution of (r+, X(r+)) 
and also for that of (t~,X(t~)). 
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Proposition 5.9 We have, for K(A) > 0, 

#J-i 



E 



-At. 



•,X(r+)€«faJ/dz = 2 \~ q / N 

q=0 



#K-1 



^ q/N 

q=0 



5®(z) forx^a, 

(5.12) 
5^\z) forx^a. 



keK 



The appearance of the successive derivatives of 8 a suggests to view the distribution of 
(t+,X(t+)) as a tempered Schwartz distribution (that is a Schwartz distribution acting 
on the space 5 of the C°°-functions exponentially decreasing together with their derivatives 
characterized by 

V<^e<S, J J ^(z)P :E {T+Gdt,X(T+)Gdz}=E :E ^(T+^(X(T+))]. 



5.2.2 Inversion with respect to A 



In order to extract the densities of (r+,X(r+)) and (r a , X{r a )) from (5.12), we search 
functions Ii q , ^ q ^ max(#7 — 1, # J — 1), such that, for §£(#/£) < 0, 

(5.13) 



/ <X> e-%(t;0rfi = A-«/V'^. 

JO 



The rhs of (5.13) seems closed to the Laplace transform of the probability density function 
of a completely asymmetric stable random variable, at least for q = 0. Nevertheless, 
because of the presence of the complex term Q\ within the rhs of (5.13), we did not find 
any precise relationship between the function I\ q and stable processes. So, we derive below 
an integral representation for Ii q . 

Invoking Bromwich formula, the function Ii q writes 

1 f* I OO 1 f~ \~ OO 

l lq (t;0 = — / \- q / N e tx+d ^d\ = — / (i\)-%e itx+e ^d\ 
' .oo 2vr J_ oq 



1 

2ivr 
1 

27 



e~2N 



r+oo in r+oo 

Jo Jo 
The substitution A i — ► X N yields 



N 
2n 





iirq j 


e 


2JV / 




Jo 




r+oo 


i-Kq 




2N 






Jo 



^N-q-l e -itX N +eie ^^dX 



I 17T 

e 2n\ together with the residues theorem provide 



and the substitutions A i- 



In particular, for q = we have, by integration by parts, 



Ni 


e 


r+oo 

mq / 

Jo 


iitq 

— e n 


r+oo 

/ X N ~ 






lo 



lio(t;0 



MA 
2vrf 



h 

IN 

Jo 



' e -tx"+e l e i $zx d x_ e -% 



r+oo 

Jo 



e -tx^+e ie -W^ dX 



.(5.14) 
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Remark 5.10 The following relation holds between all the functions ij 9 's: 



d m I, 



dp 



(t;0 = Orit q -m(t;0 forO^m^g. 



Hence, (5.12) can be rewritten as an explicit Laplace transform with respect to A: 

e - A ^ + ,X(r+)Gdzj /dz= [ 
1 J Jo 



E, 



» + 00 



e- Xt dt 



#J"1 / \ 

E E^^(*;*-«) 

q=0 \ je J / 



We are able to state the main result of this part. 

Theorem 5.11 The joint "distributional densities" of the vectors (t+,X(t+)) and 
( T ai X ( T a)) are 9iven by 

#J-i 



where 



; {r+ G dt,X(r+) G dz}/dtdz = J q (t;x- a)S$\z) forx^ a, 

q=0 
#K-1 

,{t~ G dt,X(T~) G dz}/dtdz = K q (t;x-a)Si q \z) for x > a, 

9=0 

^(*; = E ^(*! and ^(*; = E ^fcgfa £)■ 



(5.15) 



Remark 5.12 Another expression for J q (t;Z), for instance, may be written. Indeed, for 
£ < and < q < # J - 1, 



^(t;£) = 



iVi 
2^ 



— e n 



r+oo ( 

w (E 

^ Vie J 



Cj q Aj e 9: > eW ^ x \\ N ~ q ~ 1 e~ txN dX 



+oo 



Y J Cj q A j e e ^^)\ N - q - l e- txN dX 



(5.16) 



The second integral displayed in (5.16) is the conjugate of the first one. In effect, by 
introducing the symmetry a : j £ J i — ► a(j) G J such that = 0j, we can see that 



A 



n 



(9j 



An 



and 



C*(j)q = O-q l€J\ {cr(j)}) = 0~ q , I G J \ {j}) = G q (0J, I G J \ {j}) = Cj,. 

So, the sum lying within the second integral in (5.16) writes 



E ^)<A0) e°'W** = E <*A = E 
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As a result, 



N 

J g (t;0 = — ^ 

7T 



+00 



E CjqAj e djeW ^ x J X N ~ q ~ 1 e~ txN dX 



In particular, J q (t;£) is real and for q = we have, since Cjo = 1 and YljejAj = 1> 



N 

Jo(t\0 = — 9 

7T 



TTt 



+°° / in \ 

+00 / 
Vie J 



which is nothing but Pz{r+ G dt}/dt. 



5.3 Distribution of the hitting places 



We now derive the distribution of the hitting places X{t£) and X(r a ). To do this for 
X(t+) for example, we integrate (5.15) with respect to t: 



/"TOO 

; {X(r+) G <fe}/cfe = / P x {r+ G di, X(r+) G dz}/d2 
JO 

r /-+oo "I 
= £ / J q (t;x-a)dt 5<?\z) 



iV 



= > F _ 



g=0 



JO 



1 i7rg 

9 e - ^ 



J6J 



dA 



#>(*). (5.17) 



We need two lemmas for carrying out the integral lying in (5.17). 



Lemma 5.13 For any integers m, n such that 1 ^ n ^ m— 1 and any complexes a±, . . . ,a m 
and 61, . . . , b m such that K(6j) ^ and 9^ Sj=i a j^j) = /or ^ / ^ n — 1, 



+00 



91 > "a,-e-^ A 



Proof. We proceed by induction on n. 

For n = 1, because of the condition 9 [Y^jL\ a i) = ^' we can re Pl &ce ^( a m) by 
Y^f=i a j) • This g ives 



r +oo / rn \ , x . +00 [m-1 

The foregoing integral involves the elementary integral below: 



dX 

T' 



) y = 9 (log6 m - log&j) . 
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Therefore, 

c+oo 



f 

Jo 



m— 1 

^ a,- (log 6 m - logfe,-) 
j'=i 



E «i lQ g b 3 
3=1 



which proves Lemma 5.13 in the case n = 1. 



Assume now the result of the lemma valid for an integer n ^ 1. Let m be an integer 
such that m ^ n+2 and ai, . . . , a m and b±, . . . , b m be complex numbers such that $t(bj) ^ 

and 9^ Sjli a i^j) = for ^ Z ^ n. By integration by parts, we have 



E 



rfA 



~| OO 

-5>(J^)J o -;jr a (£^)^ 

Applying L'Hopital's rule n times, we see, using the condition 3^ Sj=i a i^j) = f° r 



^ Z ^ n, that 



-^r^ (EJli«i e_bjA ) J = °- Putting aj = ctjbj, we get 



+oo 



E 



dA 

A n+1 ~~ ~ra 7 



E 



cZA 
A™ ' 



We have Q^E^a^-) = 9 ( Ejli = 0for0<Un-l. Then, invoking the 



recurrence hypothesis, the intermediate integral writes 

»+oo 







E 1 ^ 



and thus 



E a i e ~ 



(n-1) 
M \ dA _(-!)«+! 



A n+1 



n! 



3 J^logfr,- 
\j=i 



which achieve the proof of Lemma 5.13. 



Lemma 5.14 We have, for ^ p < q ^ #J — 1, 



E c jq$jAj 



ifp^q-1, 
(-1)* i/p = </. 



Proof. Consider the following polynomial: 

/ \ #J-i 



E f E^iV*)' = E^ E 



e^- n c 1 -^) 

jeJ «eJ\{j} 

leJ jeJ J 



We then obtain, due to (2.11), if p ^ # J - 1, 

E E^%- 

g =0 VjrGJ 
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-x) 9 = x p 



which entails the result by identifying the coefficients of both polynomials above. ■ 
Now, we state the following remarkable result. 



Theorem 5.15 The "distributional densities" of X(t£) and X(r a ) are given by 



,{X( T +)edz}/dz = ^(-l)*VL-pLsM(z) forx^a, 



q=0 
#K-1 



(5.18) 



t {X(r-)€dz}/dz 



E (-i) 4 



[x — a) 



■5 { a q \z) forx^a. 



q=0 



It is worth that the distributions of X(t£) and X(t~) are linear combinations of the 
successive derivatives of the Dirac distribution 5 a . This noteworthy fact has already been 
observed by Nishioka [17, 18] in the case N = 4 and the author spoke of "monopoles" 
and "dipoles" respectively related to 5 a and 5' a (see also [19] for more account about 
relationships between monopoles/dipoles and different kinds of absorbed/killed pseudo- 
processes). More generally, (5.18) suggests to speak of "multipoles" related to the 

In the case of Brownian motion (N = 2), the trajectories are continuous, so X(t^) = a 
and then we classically write ¥ x {X(t^) G dz} = 5 a {dz) where S a is viewed as the Dirac 
probability measure. For N ^ 4, it emerges from (5.18) that the distributional densities 
of X(Tq) are concentrated at the point a through a sequence of successive derivatives of 
S a where 5 a is now viewed as a Schwartz distribution. Hence, we could guess in (5.18) 
a curious and unclear kind of continuity. In Subsection 5.6, we study the distribution of 
X{t^— ) which will reveal itself to coincide with that of X(r^). This will confirm this 
idea of continuity. 

Proof. Let us evaluate the integral lying in (5.17). We have, thanks to Lemma 5.14, 

e-^E ^ (% ef ) 1 = ef (l - q) E ~^ A A = if u ? - 1. 

ieJ jeJ 
Therefore, the conditions of Lemma 5.13 are fulfilled and we get 

P*{X(r+) G dz}/dz 



= E 

</=0 

= E 

q=0 



(-l) q N 
irq\ 

(-l) q N 
irq\ 



E 9 (^'^log(%ef)) 



ieJ 



x — a\ 



» ( E ] + jf » ( E 



The second sum lying within the brackets is equal, by Lemma 5.14, to (— l) q . The first 
one vanishes: indeed, by using the symmetry a : j G J i — ► a(j) G J such that 9 a ^ = Qj, 

® ( E (^M ? ) ar ^) ) = 5 ( E ^M" ar g(^) + E c n^) arg(<9j) J 
Vie J / Vie J ieJ / 

= ^ ( E s i^^l arg(0j) + ^(,)5^(J)^(,) ar S^-0)) ) " 
VieJ ieJ / 
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The terms of the second last sum are the opposite of those of the first sum since 

c a(j) q A a (j)O q a{j} = CjqAjO] and arg(0 ff(j -)) = - arg(^) 
which proves the assertion. As a result, we get (5.18). ■ 



5.4 Fourier transforms of the hitting places 

By using (5.18) and (5.11), it is easy to derive the Fourier transforms of the hitting places 
X(t+) and X(t~). 



Proposition 5.16 The Fourier transforms of X(t£) and X(r a ) are given by 



E, 



t/ia 



#J-1 

E 

q=0 
#K-1 

E 

q=0 



(x — a) 



[x — a) 



(i/J,) q for x ^ a, 



(ifi) 9 for x ^ a. 



(5.19) 



In this part, we suggest to retrieve (5.19) by letting A tend to + in (5.4). We rewrite (5.4), 
for instance for x ^ a, as 



E. 



-\T+ +ilJ,X{T+) 



n 

leJ v 
• \ #J-i 



N/T - I Z_— / -i ill 



IF* r"° 



xu(i " jA ' '" Xx ' 



^ / ?eJ 1 : 



(5.20) 



Using the elementary expansions, as A — > + , 



#J-1 



A (#</-D/^y 



^(*-a) = i. ^.VA(x-a))%o(A(# J - 1 )/ JV ), 



yields 



E 



jA? ^-^(x-a) 



?eJ 1 - 



eJVx 

ifi 



jeJ 



#J-l / r 

E E 

r=0 \g=0 



(x — aY 



em 



+ 



x (#J~l)/N^j 



= *E (EC A i) IE 

On the other hand, applying (2.11) to x = gives 



x — a 



,\{ifiy 



-\ A^ + o(A(# J - 1 )/ JV ). 



if r^#J- 2, 

' 1 (-l)^- 1 ^-^^ ifr = #J-l. 
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Therefore, 



E- 

jeJ 1 



9jAj N, 



i\x 



VjeJ / V g=o 



(x - a)« 



iA (#j-i)/jv_ 



Consequently, the limit of E : 



when combining (5.20) and (5.21) is 



(5.21) 

as A — ► + ensues. The constant arising 



(-1) 



#J-i 



vie J / 



Y[e j )(-ifJ,)* J - 1 e ifxa = (i fJ ,)* J - 1 e ifM1 . 
jeJ J 



In view of (5.19), we have proved the equality 



lim E, 



-\r++i^X{T+) 



E, 



Remark 5.17 The distribution of X(t+) may also be deduced from the joint distribu- 
tion of (t+,X(t+)) through (5.12). Indeed, by letting A tend to in (5.12) and using 
elementary expansions together with Lemma 5.14, 



jeJ 



E c ii A i E ' 
jeJ p=o 

q 



9j{x-a))P 



pi 



which, with (5.12), confirms (5.18). 



(_!),(* - y/N 



5.5 Strong Markov property for 

We roughly state a strong Markov property related to the hitting times t^. 
Theorem 5.18 For suitable Junctionals F and G, we have 



E, 



= ^[F[(X(t)) ^ t<Ta± ) E x(Ta±) [G((X(t)) m ) 



(5.22) 



E x [G{(X(t + T+)) m )] 
E x [G((X(t + r-)) m )] 



g=0 



ifx^za, 



(5.23) 



E i ^ J^ E - v 



g=0 



if x ^ a. 



Proof. We first consider the step-process X n and we use the notations of Subsection 5.1. 
On the set {r+ n = fc/2 n }, the quantities F((X n (t)) 0<t<T +J and G((X n (t + r+ n ))^ ) 
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depend respectively on X n>0 , X Hjl , ... , X„ jfc -i and X n>k , X njk+1 , ... So we can set, if r+„ 
k/2 n , 



F\ K (Xn(t)) ^ t<T + n j = Fk(X nj o,X nj i, . . . ,X nj k-i) = F H: k-i, 

G((X n (t + T^ n )) t ^o) = Gk(X n ,k,X nt k+l, ■ • •) = G n ,fc- 

Therefore, 

oo 

F {( X ni t )) ^t<T+ n ) G (( X n(t + T^ n )) t ^o) = ^2 F n,k-lGn,k\M n>k _ 1<a ^M, 



n,k } ' 

fc=l 



Taking the expectations, we get for x ^ a: 



x[F((Xn(t)) ^ <T + n ) G((X n (t + T+ n )) t>0 )_ 

oo 

= ^Ex[-P , n,fc-ll{M„ ifc _ 1 <o<Af„ ifc }IEx„ ifc (G ! n,o) 
fc=l 

= E x [F((X n (t)) Q<t<T + n ) E x{<n) [G((X n (t)) m )]\ (5.24) 

and (5.22) ensues by taking the limit of (5.24) as n tends to +00 in the sense of Defini- 
tion 3.3. 

In particular, choosing F = 1, (5.22) writes for x ^ a 

/+00 
F x {X(r+) € <fe}EjG((X(t)) t>0 )] 
-00 

which, by (5.18), immediately yields (5.23). ■ 

The argument in favor of discontinuity evoked in [12] should fail since, in view of (5.13), 
a term is missing when applying the strong Markov property. 

5.6 Just before the hitting time 

In order to lighten the notations, we simply write = r a and we introduce the jump 
A a X = X(r a ) - X(r a -). 

Proposition 5.19 The Laplace- Fourier transform of the vector (r a , X(r a — ), A a X) is re- 
lated to those of the vectors (r a ,X(r a — )) and (r a ,X(r a )) according as, for 9ft(A) > and 
(i,v G R, 



E 



e -\T a +ifJ,X(Ta-)+ivA a X 



= E 



e -\r a +inX(T a -) 



= E 



e -\r a +ifj,X(T a ) 



(5.25) 



Proof. The proof of Proposition 5.19 is similar to that of Lemma 5.1. So, we outline 
the main steps with less details. We consider only the case where r a = and x ^ a, the 
other one is quite similar. 

• Step 1 

Recall that for the step-process (X n (t))t^o, the first hitting time r+ n is the instant 
t n ,k with k such that M n ^-i ^ a and M n ^ > a, and then X{r a ^ n — ) = X n ^-i and 
X(T a!n ) = I n , fc . Set A„ jfc = X Hik - X nyk _ x . We have, for x ^ a, 

e -\Ta,n+ilJ.Xn(Ta,n-)+iv&-aXn 
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E (= -At n ,A;+iM^n,fc-l+^A n ,fc-| 
e 1 {M„ ifc _ 1 ^a<M n , fc } 



fc=l 



+ 



fc=l 



• Step 2 

We take the expectation of (5.26): 



E, 



X n 



^-\/2 n +inx+K N (iv) N /2 n 



-\/2 n +ifJiA n , k +ivA 



fc=i 



The expectation lying in the rhs of the foregoing equality can be evaluated as follows: 

Ex 



J — oo 

= / e i " 1 'P x {X nifc _i G dy, M n , fc _i ^ a} 

J — oo 

x e -^ Eo ( e ^,i 1{Xn ^ a _ j/} )E ( e ^, 1 ) -E (e^l {Xn ^ a _ y} ] 
For computing the term within brackets, we need the following quantities: 

E (e^ or ^ x "- 1 l {Xnl ^ a _ J/} ) = p V e i ^ oru >p( y l/2 n --z)dz, K (e iuX ^) = e*"^/ 2 ". 

J —oo 

With these relations at hand, we get 



E, 



-\T a , n +iflXn{Ta,n. — )+ivA a Xn 



— g-(A-Njv(w 



1 f 1 

^ fc=l ■'-oo 



x2 r 



-(\- KN (w) N )/2 



1 T V e^ z p(l/2 n ;-z)dz- r V e ivz p{l/2 n ;-z) dz 

J -co J —CO 



. (5.27) 



E 



Step 3 

We now take the limit of (5.27) as n tends to infinity: 

-\T a +i^X(T a -)+iuA a X 



poo pa 

= e ilIX + e~ xt dt e ilxy <p(X,n,u;y)F x {X(t)edy,M(t)^a} 

JO J -oo 



where we set, for y < a, 
¥>(A, ^ ; y) = lim - 

e^0+ £ 



-(\- KN (iv) N )e 



a-y 



e iflz p(£;-z)dz 



/a-y 
-oo 



e ivz p{e--z)dz . 

(5.28) 



• Step 4 

For evaluating the above function 92, we need two lemmas. 
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Lemma 5.20 For ^ p < N, we have 

' 1 

E [X(ty] = < 



for p = 0, 
for 1 ^ p ^ N - 1, 

K N Nlt forp = N. 



Proof. By differentiating k times the identity Eo(V uX( ^ = e K N^u) N t reS p ec t to u 
and next substituting u = 0, we have that 



E [X(t) k ] = H) 



9 fc 



K N (iu) N t 



u=0 



Fix a complex number It can be easily seen by induction that there exists a family 

of polynomials (Pk)keN such that, for all k £ N, 

|^(e- Ar )=P fc (u)e^. (5.29) 
In particular, we have Pq(u) = 1 and Pi(u) = Nau 1 ^^ 1 . Using the Leibniz rule, we obtain 



Pk{u) 



d k / n\ n d k ~ l 



/. . | \ u N+j-k 



Nau N - l e auN ^j 



du k V 

- Nta E (TJ(^W Fj(u) ' 

j=max(0,fc-AO V 7 V J 7 

This ascertains the aforementioned induction and gives, for u = 0, 

if 1 s: k ^ N - 1, 



Pfe(0) 



iV!aP (0) = iV!a if fe = AT. 



Choosing a = K^i t and n = in (5.29), we immediately complete the proof of Lemma 5.20. 



Lemma 5.21 For a < < (3, the following expansion holds as e — ► + : 



/ e^p(e;- 

J a 



z) dz = 1 + KN(ifi) N e + o(e). 



(5.30) 



Proof. Performing a simple change of variables and using some asymptotics of [12], we 
get 

■ [V £l/N ■ 1/N /• +0 ° • 1/N 

/ e^ z p{e--z)dz = / e lfl£/ z p(l;-z)dz = / e^ e ' z p(l; - z) dz + o(e) 



= E 



l a / £ l/N 



zP /N 



p=0 



pi 



+oo 



z p p(l;-z) dz + o{e). 



Observing that J_ 00 z p p(l; —z) dz = Eq(X(1) p ), we immediately derive from Lemma 5.20 
the expansion (5.30). ■ 

• Step 5 

Now, plugging (5.30) into (5.28), it comes 

ip(\,H,v;y) = lim - [(l - (A - K N (iv) N ) e + o(e)) (l + K N (ifi) N e + o(e)) 
- (1 + KN (iv) N e + o(e)) ] = -A + n N (ifi) N . 
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Therefore, 



E. 



-\T a +inX(T a -)+ivt\ a X 



poo _ 

- (A - K N (ifi) N ) / e~ M E x e ifxX ^,M(t) ^ a 
Jo L 



dt. 



Writing finally 



J o 



er xt E x 



e^ x(t \M{t) ^ a 

poo 

Jo 



dt 



A — KN(i/J>) N 



f 

Jo 



dt 



e- xt E, 



e^ x(t) ,M{t) > a 



dt 



dt, 



we obtain (5.25) by invoking the relationship (5.1) and by noting that the result does not 
depend on v (and hence we can choose v = 0). ■ 

Choosing (i = or v = 0, we obtain the corollary below. 
Corollary 5.22 We have, for 5R(A) > and fi, v G R, 



E a; [e- Ara+ ^ X(ra " ) l = Eje- ATa+ ^ x(ra) l and Ej e - Ara+i!/A 



E x e~ Ara 



From Corollary 5.22, we expect that A a X = X(r a ) — X(T a —) = in the following sense: 
VV G S, Eo[ip(A a X)] = 93(0). This provides a new argument in favor of continuity. 



5.7 Particular cases 
Example 5.23 Case N = 3: 

• In the case K3 = +1, densities (5.15) write 

Ee{t+ G dt, X(t+) G dz}/dt dz = J {t; x - a) 5 a {z) for x < a 

and 

P x {r~ G (it, Af (t~) G dz}/dtdz = /Co(i; x — a) 5 a {z) + /Ci(i; x — a) S' a (z) for x ^ a. 
Here, we have d\\ = 62 = 9\, d<i\ = 61 = #2 and 



Jo(*;0 



irt 

7rf 



ijr /-+00 

ef / 6*0^0 e' 
■/ 

Jo 

+°° • />/3 



d\ 



d\ 



e2 



Ux-tx 3 



7T 



sin ( — £A + - ) d\; 



— y 


277 

e 3 


/ 

Jo 


vrv^t 


e 3 






+00 



B X B X e e ^^ x + 2 B 2 e ^ e¥ « A ) e- tx3 dX 
p+00 

f l ( 



e * — e 



7T 



e 3«* cos — £ A 



e" tA3 dA; 
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/Ci(*;0 = — » 

7T 



e s 



+00 



e 3 



+00 



e -€A_ e e-*^ Ae -U3 dA 



/■ +0 ° Tl fA i fA /V3 vr 



Ae""*^. 



In the case K3 = —1, densities (5.15) write 

P :r {r+ G dt,X(r+) G dz}/dtdz = J {t;x- a) 5 a {z) + Ji{t;x - a)5' a {z) for x s£ a 
and 

P x {r~ G (it, X(t~) G dz] /dtdz = JCo(t; x — a) 8 a {z) for x ^ a. 
In this case, we have cqi = 62 = do, c 2 i = Oq = 62 and 



'0 



e-* A3 dA 



A e~ dA 



e 2 



Let us point out that the functions J7b, Ji, K.$ and /Ci may be expressed by means 
of Airy functions. 



Example 5.24 Case N = 4: formulas (5.15) read here 

P x {r+ G dt,X(r+) G dz}/dtdz = J {t;x- a) S a {z) + Ji(t; x - a)5' a {z) for x < a 

and 

P x {r~ G (it, X(r~) G dz] /dt dz = /Co(t; x — a) <5 a (z) + /Ci(t; x — a) ^(2) for x ^ a. 
We have c 2 i = 6*3 = #2, c 3 i = 9 2 = 63, doi = 61 = 6 = -9 2 , dn = 9q = 9 1 = -9 3 and 



7Tt 



e 4 



(^ 2 A 2 e 9 ^ « A + #3^3 e 9 ^ « A ) e~' A4 dA 



vry^t 



f / (^- e ^) e -^ 4 dA 



2nt Jo 

^ r+00 



e« A - cos(£A) + sin(£A) 



e^dA 



dA; 



Mt-o = — 9 

7T 



e * 



2^/2 



7T 



+00 



+00 



jf " (c 2 iA 2 e 6 ^ « A + c 3 i^3 e° 3e * « A ) A 2 e^ 4 dA 

1 [ +0O (e^-e^ x )X 2 e- tx4 d\ 
Jo 

- V2 cos A - ^) 



-/ 

7T Jo 



cos(£A) + sin(£A) - e ?A 
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A 2 e" tA4 dX 



X 2 e~ tx dX 



and similarly 



t r+oc 



+00 



cos(£A) + sin(£A) - e~ 5A 
< x - cos(£A) + sin(£A) 



' tx4 d\; 



X 2 e- txi d\. 



We retrieve formulas (8.17) and (8.18) of [18]. 



5.8 Boundary value problem 

We end up this work by exhibiting a boundary value problem satisfied by the Laplace- 



Fourier transform U(x) = E, 



-\T++inX(r+) 



, x £ (—00, a). 



Proposition 5.25 The function U satisfies the differential equation 

V x U{x) = \U(x) for x £ (— 00, a) 
together with the conditions 

C/W( a -) = {i^) l e ^ a for ^ I ^ #J - 1. 



(5.31) 
(5.32) 



Proof. The differential equation (5.31) is readily obtained by differentiating (5.4) with 
respect to x. Let us derive the boundary conditions (5.32): by (5.4), 



e\A 3 



By (2.11) we see that 



7a ^ 



(<M)' 



which proves Condition (5.32). 



J4t_a. 



aiw(i-^) 



We also refer the reader to [19] for a very detailed account on PDE's with various 
boundary conditions and their connections with different kinds of absorbed/killed pseudo- 
processes. 
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